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We investigate the use of perturbation theory in finite sized frustrated spin systems by caicuiating 
the effect of quantum fluctuations on coherent states derived from the ciassicai ground state. We 
first caicuiate the ground and first excited state wavefunctions as a function of appiied fieid for a 12- 
site system and compare with the results of exact diagonalization. We then apply the technique to 
a 20-site system with the same three fold site coordination as the 12-site system. Frustration results 
in asymptotically convergent series for both systems which are summed with Pade approximants. 

We find that at zero magnetic field the different connectivity of the two systems leads to a triplet 
first excited state in the 12-site system and a singlet first excited state in the 20-site system, while the 
ground state is a singlet for both. We also show how the analytic structure of the Pade approximants 
at |A| ~ 1 evolves in the complex A plane at the values of the applied field where the ground state 
switches between spin sectors and how this is connected with the non-trivial dependence of the 
< S z > number on the strength of quantum fluctuations. We discuss the origin of this difference in 
the energy spectra and in the analytic structures. We also characterize the ground and first excited 
states according to the values of the various spin correlation functions. 

PACS numbers: 75.10.Jm Quantized Spin Models, 75.50.Ee Antiferromagnetics, 75.50.Xx Molecular 
Magnets 



I. INTRODUCTION 



The antiferromagnetic spin Heisenberg model has been the object of intense investigation through the years. Re- 
cently it has attracted enormous interest in the study of strongly correlated electron systems, which include the oxide 
superconductors Q and low dimensional spin systems. In the limit where the on-site Coulomb repulsion is very strong, 
it is equivalent to the one band Hubbard model for half filling. The inclusion of competing interactions in it has led 
to novel quantum phases, making it appropriate for the study of quantum criticality 0. 

The solution of the model was calculated by Bcthc in one dimension for nearest neighbor interactions jj), but a 
solution in analytic form is lacking for two or three dimensions, except for special cases [Q. Approximation and 
numerical techniques that have been used include diagonalization of small clusters |||| , Monte Carlo techniques , 
cluster expansions || , spin wave expansions [p|~|ri|| and the density matrix renormalization group JO] . The methods 
that consider the full Hilbert space of the problem are limited by the size of the system, since the number of states 
is exponentially dependent on it. In Monte Carlo calculations the sign problem leads to loss of statistical accuracy, 
especially for frustrated systems 0. 

An alternative approach to these techniques is the direct application of perturbation theory in which corrections to 
the classical treatment are calculated order by order in the strength of residual interactions, the effect of fluctuations 
left out in the mean field approximation. Although this approach can also be limited in the size of the systems 
which can be investigated due to the dimensionality of the Hilbert space, it is very different in that it provides 
analytic information on the effects of corrections to classical approximations and so complements the other approaches 
mentioned above. We use this approach to calculate the ground state wavefunction of two systems which have 
qualitatively different classical ground states, as we will explain below. We demonstrate that the differences survive 
in the exact spin i ground states through the difference in the analytic structures of the Pade approximants derived 
from the perturbation expansions. 

We consider the Heisenberg model on closed two dimensional spin i systems which have three nearest neighbors. 
The antiferromagnetic interaction between the spins leads to frustration at the classical level. One such system is the 
60 site cluster whose relative positions are the same as those of carbon atoms in Cqq (from now on we will refer to 
n site systems with the spatial symmetry of the fullerenes as C n ). The 60 site system consists of 20 hexagons and 
12 pentagons. Assuming a tight binding model for its electronic properties where there is one orbital per site and a 
strong Coulomb repulsion for double occupancy, the Heisenberg model on it gives an effective low energy description 
at half filling. This is a first approximation to the problem, since Cqo is estimated to be in the intermediate and not 
in the strong coupling U regime (Of . The hopping matrix elements between sites on the same pentagon can be taken 



to be different from the ones between sites on adjacent pentagons. This leads to two positive exchange constants J\ 
and J2 and the Hamiltonian is: 
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H = ,hj2 Si-Sj + JzJ^ & ■ 4 • (!) 

<i,j> <i,j> 

J\ refers to bonds between the same pentagon and J2 to non-pentagon bonds, while <> stands for nearest neighbor 
interactions. 

The classical ground state of the Heisenberg Hamiltonian has been determined on this and analogous C n systems 
fl3|| . The magnetic properties of these ground states show an unexpected dependence on n. In some the magnetization 
is discontinuous in an applied field whereas in others it is the susceptibility which is discontinuous. It is of interest 
to determine whether this dependence is an artifact of the classical approximation or whether it is present in S = h 
solutions. We use the classical C n ground states to define an Ising Hamiltonian where the quantization axis at each 
site is determined by the direction of the spin vector in the classical ground state. Coherent states are defined along 
these axes and constitute a mean field approximation for the quantum ground state. In this way, each site has a 
local axis associated with it, in contrast to the basis where all spins are defined in the same coordinate system in 
spin space. From now on we call the former 'local basis', and the latter 'global basis'. The quantum fluctuations 
are then built around the classical directions. They are the terms added to the original mean field approximation, 
and they are multiplied with a parameter A which varies from to 1. When A — 1, the full isotropic Heisenberg 
Hamiltonian is recovered. Thus the solution is generated as a series expansion in the perturbation parameter A with 
application of perturbation theory, and the expansions of the ground state energy and wavefunction in the local axes 
basis are known. With this approach we can study the evolution of the system away from the classical ground state 
and towards the full quantum limit A = 1. 

However, the Hilbert space for the 60 site system is huge, consisting of 2 60 = 1.15 x 10 18 states, so a perturbation 
treatment in the whole Hilbert space can only give a few orders. This is because the number of states involved 
in the calculation rapidly increases as the order and frustration increase ]l4| . In addition, since the spin axes are 
directed along the classical solution's directions, the total spin in the global basis is not a good quantum number of 
the Hamiltonian (when A < 1), so a reduction of the number of states by focusing on a particular S z value is not 
possible. Another way of gaining insight into the problem is to consider similar smaller systems, belonging to the 
same family. In all these systems there are rih = § — 10 number of hexagons and 12 pentagons. The smallest member 
of the group is the 20 site system. Again, we consider S — | spins sitting at their vertices. 

The frustration of the Hamiltonians studied leaves its signature in the various series expansions generated by 
perturbation theory, producing functions with non-analytic structure in the complex A plane. The presence of branch 
cuts limits absolute convergence within a circle centered at the origin with a radius of convergence smaller than 1 and 
the series are only asymptotically convergent in the full isotropic limit where A = 1. Therefore, we employ analytic 
continuation with the use of Pade approximants. The structure of the functions in the A plane depends strongly on the 
form of the perturbing Hamiltonian. We investigate the signature of the changes in the complex plane structure of the 
related functions as the perturbing Hamiltonian is varied. The generalization to complex variables has been proven 
useful in the study of phase transitions in the two dimensional Ising model in temperature T in a complex magnetic 
field jl{|. The systems studied here are closed and do not possess a thermodynamic limit. However, knowledge of the 
structure in the complex coupling constant plane provides information about the functions studied p6| . It can also 
provide information for the evolution of instabilities. Here the only possible transitions are between < S z > sectors 
with increasing magnetic field, where < S z > is the expectation value of the z component of the spin in the global 
basis. 

The classical ground state at zero magnetic field is doubly degenerate, since a flipp ing of all the spins does not 



change the energy. Consequently degenerate perturbation theory has to be applied |17|. This is done via a similarity 
transformation, and an effective 2x2 matrix is generated. This matrix provides information on the ground and first 
excited state, including the evolution of < S z > as a magnetic field is varied. The structure of this matrix in the 
complex A plane can be correlated with the transitions of the < S z > number between different sectors as a function 
of the magnetic field. 

The plan of this paper is as follows: in section II the method for the solution of the problem using perturbation 
theory is described. In section III this is applied to a twelve site system, C12, where the results are compared with exact 
diagonalization and found to be in complete agreement. We also discuss how the analytic structure of perturbation 
theory is reflected in the magnetic properties. In section IV perturbation theory is applied to C20 and the results are 
tested by recovering expectation values for S 2 and S z very close to integer values. Here we contrast the magnetic 
field dependence of the analytic structure with the results found in the classical approximation. For both C12 and 
C20 ground and first excited state wavefunctions are calculated in applied magnetic fields and the dependence of their 



magnetic properties on the strength of quantum fluctuations is determined. In the case of C20, the ground and excited 
states at zero magnetic field are singlets. This result is in agreement with similar ones for strongly frustrated magnetic 
systems such as the Kagome lattice and a one dimensional analogue of the pyrochlore lattice Eq] . The non-magnetic 
nature of the excitation is attributed to the frustration and the connectivity of the system. 



II. METHOD 

The starting point in the calculation is the classical ground state. The Hilbert space is spanned by spin ^ spinors 
determined by the classical solution. Local Zi axes are defined along the classical spin directions, and spin states are 
defined at each site such that the expectation value of the component of the spin along the axis equals its classical 
value: 

where on, fii are spinor coefficients and \Si \ = 1. The coherent states are products of spin states along these axes: 
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where N is the number of spins in the system. Here we limit ourselves to these 2 N states out of the overcompletc 
basis of the coherent states. These are eigenstates of the unperturbed Hamiltonian Ho and constitute an orthonormal 
basis. The classical ground states are the ones where all the spins are either 'up' or 'down' so that Hq has the form 
of an Ising Hamiltonian with respect to the local quantization axes defined by the classical result. The fluctuations 
around the local axes are raising and lowering spin operators, defined along the local a;, and j/i axes. These are 
included in the perturbing part Hi which we scale with a parameter A. Each a;, axis is defined in the azimuthal plane 
of the corresponding axis and perpendicular to it, and then the y$ axis is constructed with the right hand rule. Due 
to the absence of a global z axis, the Hamiltonian takes a complicated form. The solution is generated in perturbation 
theory as a power series in A. Therefore the Hamiltonian has the form: 

H = H a + \H 1 , (4) 

where Hq is the classical part and Hi the quantum fluctuations. 

There are two well known issues which complicate the application of perturbation theory. The first is the double 
degeneracy of the coherent states generated from the classical ground state in the absence of an applied magnetic 
field. This requires the use of degenerate perturbation theory. To overcome this problem, an effective Hamiltonian is 
constructed for the degenerate ground states, via a similarity transformation fTfl |. The Hamiltonian is transformed 
to a block diagonal form and its elements, as well as the coefficients of the Hilbert space states, which contribute to 
the perturbed wavefunctions, are expanded as power series in A. Then recurrence relations can be written down for 
the effective Hamiltonian and the wavefunctions: 

H e k ff (m,£) =< ml-Hil*^! > , (5) 
< ^ >= Eo~^ {< n|i/l| *^i >-EE H e k f _{,(£',£) < n\¥ k 9 >) . (6) 
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Here k is the order of perturbation, |m > a degenerate ground state of Ho, \n > a state in the Hilbert space different 
from \m >, E n its energy at the classical level and Eo the classical ground state energy. L is the dimensionality of 
the degenerate subspace (here L = 2) and i runs from I to L. The result of this calculation is an L x L (here 2x2) 
matrix whose eigenenergies are the ground state and the first excited state. This method can be applied for zero or 
non-zero magnetic field, where all the magnetic field terms in the Hamiltonian are included in the perturbation. If 
the magnetic field is included in Ho non-degenerate perturbation theory can be applied in this case. 

The second issue complicating the problem is that the perturbative part of the Hamiltonian generates series expan- 
sions whose radius of convergence does not extend to the isotropic limit A = I but rather is limited to a circle with 
a radius smaller than 1. Therefore, we analytically continue outside the radius of convergence with the use of Pade 
approximants. These are described in detail in the literature pof . The algorithms used for their calculation here are 
the determinant algorithm and the Viscovatov algorithm plj |. We now apply the method to two systems, C12 and 
Cio- 



III. Cl2 



A smaller system with a geometry similar to the ones considered above is the truncated tetrahedron C12 p2f , shown 
in figure 1. This is made up of four triangles and three hexagons, and looks like a closed triangular lattice. Every site 
is three fold coordinated, and bonds between sites on the same triangle (tr.) are called J\, while different triangle 
bonds (n.tr.) are named J2. The exact solution of the Heisenberg Hamiltonian has been found for this system and 
is used to check the results from perturbation theory here. At the classical level it was found to have a jump in the 
susceptibility at a critical value of an external magnetic field and so has similar properties to the systems discussed in 
reference [13 ]. As mentioned above the classical solution determines local axes of quantization Zi and determines Hq. 
Spins belonging to the same triangle are coplanar at the classical level and each local Xi axis is defined in this plane. 
With application of the right hand rule the corresponding yi axis is defined. With a = ^ the Hamiltonian becomes: 

^-=H Q + XH 1 , (7) 

where 

tr. 1 n.tr. 
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where cf> — tan _1 (2\/2). The coefficients reflect the dependence of the Hamiltonian on the classical spin directions. 
The Hamiltonian is complicated, since it includes any possible combination of raising, lowering and S? operators, 
and has complex coefficients. There is no choice of the local coordinate systems that would make all the coefficients 
real, once the local z axes are fixed along the classical solution's directions. The local Xi and yi can also be defined 
without reference to the specific form of the classical solution, as was stated in section II, but this results in a more 
complicated expression for the Hamiltonian. 



A. Ground State Energy and Wavefunction 



The ground state energy and the wavefunction coefficients are functions of A and a. The a — case corresponds to 
four isolated triangles, while the a — > 00 case corresponds to spins forming singlets (dimers) via the Ji bond. There 
is no further frustration when assembling the tetrahedron from the individual triangles, since this costs nothing in 
energy for the classical spins. Consequently for this system the classical ground state is independent of a and the 
quantum fluctuations select a unique ground state when A ^ 0. This is reminiscent of the order-disorder transition 
induced by quantum fluctuations for frustrated systems [ p3[ . The effective Hamiltonian for the two degenerate ground 
states has the following form: 

■ (10 » 

where An(\), Bn(\), Cjv(A) and Dn(\) are polynomials in A of Nth order, the order of perturbation. A can assume 
complex values, since we are also interested in the structure of the functions in the complex plane. The coefficients 
of series Cjv(A) and -Djv(A) are complex conjugates of the coefficients of Bn(A) and ^at(A) respectively. Therefore, 
in the case of physical interest where A is real, CV(A) and Dn(X) are the complex conjugates of Bn(X) and Ajv(A) 
respectively. Diagonalization of the matrix in the latter case shows that its eigenvalues are real, as expected. The two 
classical ground states, |0) and |0), evolve in the following manner as functions of A: 

N d-L 
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A 



and similarly for l^g)^, where Fjy-(A) are polynomials of Nth order in A with complex coefficients, \i > is a Hilbert 
space vector different from |0) and |0) and d is the dimensionality of the Hilbert space. L is again the dimensionality 
of the degenerate subspace (here L = 2). 

1. Analytic Continuation 

The next step is to analytically continue the polynomials ^4jv(A), Bn(X), Cn(X), £>at(A) and F^(X). Then the 2x2 
matrix is diagonalized and the energies as well as the wavefunctions are known. This gives the ground and the first 
excited state at each value of A. For C\i the dimensionality of the Hilbert space is d = 2 12 = 4096. 

The ground state energy for different a's is given by the converged value of Efj(X), which is the lowest eigenvalue 
for the TVth order approximant. The criterion for convergence is that at order N mean square fluctuations should be 
1% of the mean value for the 7 approximants placed around the Nth order (i.e.,iV — 3, N, N + 3). This is easily 
satisfied for small values of N when A < 0.5 for a = 1 , but the mean square fluctuations increase for larger values 
of a. The dependence on iV for different values of a can be understood by looking at the structure of Ajv(A) in the 
complex coupling constant plane (figures 2 and 3) which is discussed below. The analytically continued form for the 
polynomials is given as the ratio of two polynomials defined in the complex A plane so that, for example: 

A N (X)^A^-(X) = ^^- , (12) 
Qn(X) 

where a.c. stands for analytically continued. For A on the real axis D^ c - = (A^ c -)* and C^ c ' = (B^ c ')* guaranteeing 

real eigenvalues for real A. The structure is revealed by looking at zeros and poles of A^ C (X) for complex A. These are 

the roots of Pjv(A) and Qn{X). Frustration leads to branch cuts in the A plane which are given by lines of mixed zeros 

and poles of A^ C -(X). This is how the single valued A^ C -(X) tries to reproduce a multi- valued function associated with 

a cut in the complex plane. This structure is shown for a = 1 in figure 2 and for a = 2 in figure 3, for N = 220. The 

radius of convergence A c (circles) shrinks as a increases and for A > A c the perturbation expansion is asymptotically 

convergent. This indicates that more orders are needed for Pade approximants to converge as the ratio j- increases. 

One explanation for this is that, as a increases, the classical ground state is not as good starting point as for small A, 

since the spins tend to form singlets via the J 2 bond in the quantum limit. Thus it gets harder to reach the quantum 

state from the ordered classical ground state, the former made up of independent singlets on the non-triangle bonds 

at the limit — * 00. 
Ji 

2. Numerical Precision 

As we go to higher orders in perturbation theory the number of calculations goes up, increasing the possibility for 
significant propagation of numerical error p3|. To extend the calculation to higher orders, the package MPFUN was 
used p5| , which allows arithmetic to very high precision, limited only by machine specifications. The perturbation 
and analytic continuation were done by using typically 94 digits precision in MPFUN, except in some cases where 
194 digits were used. 

To reduce memory requirements and execution time, perturbation theory was first run in double precision. After 
the states with equal coefficients or coefficients differing only in sign (in real and/or imaginary part) due to symmetry 
were identified, the program was run with MPFUN, taking advantage of these symmetries. Thus the scale of the 
calculation was significantly reduced. The time required to get 250 orders in perturbation theory was approximately 
thirty five minutes, when eight processors were used in parallel on a SGI machine. These were IP27 processors with 
a frequency of 250 MHz. We typically used 250 orders to get convergence for fields smaller than J\. 

A second criterion of convergence comes from the total spin component < S z >, which is a good quantum number 
for the Hamiltonian in the isotropic case where A = 1. If the analytically continued wavefunction coefficients converge, 
then the calculation of < S z > should yield an integer in an applied field. Any deviation from an integer value indicates 
lack of convergence, which could be due to insufficient orders of perturbation theory used or propagation of numerical 
error, and indicates that higher numerical accuracy is needed. This criterion is more stringent since now all the 
coefficients F]^(X) have to converge but once convergence is achieved any expectation value can be calculated, since 
the knowledge of the wavefunction fully solves the problem. Since all the states in the Hilbert space determined by 
the local z axes contribute to the ground and first excited state coefficient functions, the calculated wavefunction is 
essentially exact. 



3. Degeneracy 



In the absence of an applied field states with different < S z > values are degenerate and the excited state, a triplet 
in C12, is a linear combination of S z = 0, ±1 states. However it is found that starting in the S z = sector as we do 
the perturbed state stays in this sector for all A for both ground and excited state at h=0. Although this degeneracy 
is lifted by the applied field, there is an additional degeneracy due to the geometrical symmetries of Cyi- The exact 
diagonalization shows that there are three degenerate triplet states at the first excitation energy so that even in 
an applied field the excited state with S z = 1 is triply degenerate. Choosing a particular classical ground state to 
generate the starting states |0) and |0) picks out a linear combination of these degenerate states. 

4- Correlation Functions 

After the calculation of the ground state wavefunction, its correlation functions can be directly evaluated. In Cn 
there are five kinds of qualitatively different correlation functions, two of which refer to nearest neighbors. In figures 
4 and 5 these are plotted as a function of the strength of the quantum fluctuations A for the case a = 1 for the ground 
and excited states respectively. The magnitude of the correlation function < Si ■ Sg > is smaller than 0.02 and is not 
plotted. For A = 1 the solution of the full isotropic case agrees with the one found from exact diagonalization p2| . 
The nearest neighbor correlation functions are < Si ■ S2 >— —0.125 and < Si ■ S4 >= —0.250 at the classical level 
where A = 0. Spins 1 and 4 are antiparallel, while < Si ■ S2 > is one half of < Si ■ §4 >. 

In the ground state (figure 4) < Si ■ S2 > increases in magnitude with A and at A = 1 equals —0.183, being roughly 
50% bigger than its classical value. < Si • S4 > increases in magnitude as the quantum fluctuations become stronger 
and reaches its maximum just before A = 0.8. Then it doesn't change significantly and its magnitude at A = 1 is more 
than two times larger than its classical one, being equal to —0.586. This is because as the perturbation is turned 
on adjacent spins belonging to different triangles want to create singlet bonds due to their J 2 interaction. The value 
of < Si ■ S4 > at A = 1 is close to the value for a singlet state between two spins, which is —0.750. The other two 
correlation functions have a non-trivial dependence on A, but their values in the isotropic case do not significantly 
differ from their classical ones. 

The excited state correlation functions of figure 5 do not differ significantly from the ones of figure 4. However, 
we observe a smaller value for < Si ■ S5 > compared with the ground state value, while < Si ■ S4 > has a smaller 
magnitude at A = 1 again compared to the ground state value indicating a more triplet character for this bond. 

In figures 6 and 7 we plot the correlation functions for a = 2. We observe that for both the ground and the excited 
state all the correlations decrease in magnitude, except the ones between neighboring spins connected via the J2 bond. 
The value of < Si ■ S4 > in the ground state is now —0.697, approaching the singlet value —0.750 even closer. 

B. Magnetic Field 

We next introduce a magnetic field in the problem and study the ground and excited states as a function of A. 
There are two possible approaches to this calculation. The first is to perturb around the zero magnetic field classical 
ground state. The second is to calculate the classical ground state in the presence of a magnetic field, and then apply 
perturbation theory. Because the field breaks the time reversal symmetry, the classical ground state is non-degenerate, 
and in the second case non-degenerate perturbation theory is used. The results of the two methods should agree at 
the isotropic limit, A = 1. However, their dependence on A is different. 

1. Degenerate case. Field independent classical ground state. 
The Hamiltonian for the first method has the form: 

^- = H + X(Hi + H 2 ) , (13) 
J\ 

where Hq and Hi were defined before and H2 is the part that relates to the magnetic field. Specifically: 

JV , N 

H 2 = -hY / cosOiSI + - J2 sin9i(S+ + Sr) . (14) 



h is the strength of the magnetic field in units of J\ and 0* the classical solution's angles with the global z axis. The 
direction of the magnetic field is taken along the global z axis, so it is perpendicular to the plane of one of the four 
triangles of the system. The spins belonging to any triangle are coplanar in the classical case (in spin space), and 
they can be chosen to lie in the physical plane of the triangle. 



2. <S Z > 



For A = 1, the ground state lies in the S z = sector when there is no magnetic field. < S z > is defined along the 
global z axis, and it commutes with the Hamiltonian (13) at A = 1. As the field is turned on the energy of the ground 
state, which is a singlet, won't change at the isotropic limit. However, the triplet excited state has S z = 1 there, 
and its energy decreases linearly with the magnetic field due to the Zeeman term. At a critical value of the field the 
triplet state energy becomes equal to the one of the singlet state, so the triplet state becomes the ground state. As 
the magnetic field increases further the ground state moves towards spin sectors with higher value of S z , until the 
magnetization saturates. The results for the energies found from perturbation theory and analytic continuation are 
found in figure 8, and they reproduce the exact values (A = 1). By extrapolation of the straight lines which give the 
ground and the first excited energy back to zero field, we recover the energies of excited states in the zero field case. 
The diagonalization of the effective Hamiltonian matrix gives the following eigenenergies: 



A%c-(\) + D^-(\) 



B a N 



Aa.c 



(A) 



D a -' 

N 



'(A) 



(15) 



where A^ C -(X), £% C (A), C^ C (A) and Z?^ C (A) are the analytically continued values of the corresponding polynomials 
of Nth order. At the critical fields at which the two states exchange 'roles' as the ground and the first excited state, 
the square root goes to zero at A = 1. This occurs when Bfj c - and the imaginary part of A^ c ' go to zero. In this way, 
we start from two degenerate ground states and we end up with two degenerate states. 

The strict criterion for successful convergence of the wavefunctions is the calculation of the < S z > number. This 
should be an integer at A = 1. We calculate its evolution with A in the global spin basis, where the quantization axis 
is the same for all spins. The expression for its expectation value < S z > is: 
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where G is given by 
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The plus sign corresponds to the ground state and the minus sign to the excited one. The magnitude of this coefficient 
is 1 for any real A, since in this case D^ c - — (A^ ')* and C^ c - — (B^ c -)* . The accurate calculation of the wavefunctions 
at the critical fields h c involves taking the limit liniA—i G(A). In this case the magnitude of B^ C (X) as well as the 
imaginary part of A^ c (\) go to zero, thus both the numerator and denominator vanish. The success of the analytic 
continuation is such that the calculation of A^ C -(X) and B^ c (\) is so accurate after analytic continuation that the 
result lim>_^i G(A) = ±1 is recovered and this is reflected in the calculated < S z >. The calculation was done for 
various fields for a — 1, and the results are shown in figures 9 and 10. There it is shown that indeed < S z > assumes 
integer values when A = 1, confirming the success of analytic continuation. For example, when h = 0.7, at A = 1 
< S z >= 0.99999999998 for the ground state and < S z >= 1.0 x 10~ n for the excited state. 

As seen in figure 9, the quantum fluctuations raise the value of < S z > for small values of the magnetic field, but 
eventually it goes to zero at the isotropic limit. However, just above the transition to the S z = 1 sector, which occurs 
at h c \ = 0.6878 in agreement with the exact answer, < S z > has a non-trivial behavior for intermediate A's before 
assuming the value equal to 1. Thus for h = 0.688 there is a rapid change in the magnetization at a value of A close 
to 1. There is competition between the and 1 spin sector for the ground state, and finally the quantum fluctuations 
lead to the latter. The magnetic field terms of the perturbing Hamiltonian favor a magnetized ground state, while 
the rest favor zero spin. The sudden change in < 5* z > as a function of A moves closer to A = 1 as h — > h c from 
below. If we go further away from this critical field, the jump is pushed towards smaller values of A, and eventually 
it vanishes. Similar effects are observed in figure 10, which includes the corresponding graphs for the excited state. 
In this case, the excited state has S z = 1 below and S z — above the transition. The conclusion is that there is a 
'window' around the critical field where the A dependence of < S z > is very strong. 



As the magnetic field gets bigger, more terms are needed to analytically continue the wavefunction so that the spin 
number assumes the proper integer values. There is also a need for greater numerical precision, due to the increased 
number of calculations which tend to propagate the numerical errors pifl , and so 194 digits of MPFUN accuracy were 
used. The transition between the S z = 1 and S z = 2 spin sectors for the ground state takes place at h c % — 0.9869. 
For a magnetic field equal to 1.01, 501 orders were generated to get convergence. This requirement for more orders in 
perturbation theory and more accuracy makes the calculation of the factor G of equation (16) harder as the second 
critical field is approached. The behavior of < S z > as a function of A is shown in figure 11, where the < S z > value 
of the ground state remains constant over a range of A, and eventually "jumps" to the final value < S z >= 2. On the 
other hand, the excited state < S z > approaches 2 close to A ~ 1 only to settle at 1 when A = 1. 

Non-monotonic behavior of the quantum number S z near a transition is observed for other values of J2 as well. 
These are shown in figures 12 and 13, for J 2 values equal to 0.8 and 1.8. For the second case the higher value of J2 
stabilizes the spin value for A's close to 1. It was observed that convergence was harder to get as J 2 and/or h were 
increased, for the reasons already mentioned. 

3. Correlation Functions 

Although the calculated eigenstates at A = 1 are independent of magnetic field they do depend on the field for 
intermediate values of A. This can be seen clearly in < S z > as a function of A in figures 9 and 10. As pointed 
out previously the choice of the starting classical ground state picks out a particular linear combination of the three 
degenerate (S = 1, S z — 1) states. Different linear combinations lead to different values of the < Si ■ Sj >. However 
when adding the values of nearest neighbor < Si ■ Sj > for different choices of starting classical ground states the 
values of the energies, calculated directly in equation (15), are recovered. 

A further check on the two point correlation functions calculated for the wavefunctions is < S 2 >= J2i j < Si- Sj >. 

Calculating < S 2 > in the ground and excited state at h = 0.7 we find 2.00188 and —0.0028361 consistent with 
< S 2 >— S(S + 1) for 5 = 1 and S — 0. This measurement tests the accuracy of all two point correlation functions. 
The ones for nearest neighbors for h = 0.7 are plotted in figures 14 and 15 for the ground and excited states. Comparing 
these with the A dependence of the states in figures 6 and 7, the only difference is a non-monotonic dependence at 
A ~ 0.5 especially for the nearest neighbor correlation functions < Si ■ S2 > and < Si ■ S4, >. 

4- Analytic Structure 

The motion of the zero of the square root appearing in equations (15) and (17) in the complex A plane can be traced 
out with the help of the Pade approximants, which are ratios of two polynomials. Since the analytic continuation was 
quite successful in the vicinity of A = 1, we can look for the structure of approximants in this area. To this end, we 
consider the function 

flr(A)^(A)+ ^^ (A) ; " 1A) j , (18) 

which is inside the square root in equations (15) and (17), and we plot the zeros of the numerator and the denominator 
of the approximants as a function of the magnetic field. The approximants were generated with 220 orders of 
perturbation theory (the order of the numerator and the denominator of the Pade approximant is one half of this 
number). In figure 16 the complex plane structure is plotted for a magnetic field h=0.688. The Pade approximant 
is seen to reproduce the branch cuts of the function in the complex plane. There are also zeros and poles in the 
complex plane not falling on any of the branch cuts. These are present only in approximants of specific order and 
they are artifacts of the analytic continuation. It can be seen how the zeros of equation (18) appear in the complex A 
plane close to A = 1. The coefficients of the expansion of (18) are real, so the zeros and the poles appear in complex 
conjugate pairs. 

Another feature of the picture is that it is asymmetric with respect to the imaginary A axis. This asymmetry gets 
more and more pronounced with increasing field. This is because negative A's correspond to different ferromagnetic 
couplings between neighboring sites, and these interactions favor aligned spins, a state similar to the one favored by 
the magnetic field. In this case it is more difficult to get convergence because of the competition of the classical 
antifcrromagnetic interactions with the quantum ferromagnetic ones and the magnetic field. Consequently the radius 
of convergence is smaller. 



o 



In figure 17 the roots of the approximants are plotted for various magnetic fields close to h c \ = 0.6878 in the 
vicinity of A = 1. The zeros of the square root approach A = 1 in conjugate pairs and they finally hit the axis at the 
transition. As the strength of the field is further increased, the roots move away from A = 1 in the opposite direction. 
We see this behavior in a small range of the fields around the transition with roots around A = 1, which is reflected 
in the rapid non-monotonic variation of < S z > around the critical field. Since the square root of function (18) has 
to eventually be taken, its zeros will turn into branch points in the complex A plane. 



5. N on- degenerate case. Field dependent classical ground state. 

The alternative approach to the problem is to calculate the classical ground state in the presence of the magnetic 
field, and then introduce the remaining terms as perturbation. In this case the Hamiltonian is: 

^=H +H' + X(H 1 +H' 2 ) , (19) 

where: 

N 

H' = -hJ2 cos9iS i ' ( 20 ) 
»=i 



ff; = !x>nW+sn ■ (21) 

i=l 

and Hq and Hi have already been defined in equations (8) and (9). Since the magnetic field breaks the degeneracy 
of the classical ground state Rayleigh-Schrodinger perturbation theory can be applied directly for the wavefunctions. 
The results after the analytic continuation are shown in figure 18. This method can be compared with the one of 
the previous section only at A = 1, where in both cases the isotropic Heisenberg Hamiltonian is recovered. Starting 
from the field dependent classical ground state is not any better as far as computer memory and execution time is 
concerned, since symmetries reduce memory requirements at the same level in both cases. However, the first method 
also gives the ground and first excited state very accurately, consequently the approach of the transitions between S z 
sectors as a function of the applied magnetic field is easily seen. On the other hand, the second method works better 
for higher fields, since now the classical ground state is related to the strength of the magnetic field, thus it has a 
non-zero value for < S z >. Various fields up to h = 1.3 where employed in the calculation. The results for the energy 
agree with the exact diagonalization answer . 

Looking at figure 18, we see the dependence of < S z > on A is not monotonic for the various magnetic fields. This 
can be attributed to the fact that the starting point is now the magnetic field dependent classical ground state. The 
terms in the Hamiltonian related to the magnetic field tend to increase the spin, while the rest favor a zero spin 
quantum state, thus there is competition between the two at the classical and quantum level. In the degenerate case 
the field dependent terms entered only in the perturbing part of the Hamiltonian, making the perturbation stronger 
and dominant in the determination of S z for higher A's. 

The structure of the analytic continuation of the energy was also studied in the non-degenerate case, and it was 
consistent with the structure found from the degenerate perturbation. For a magnetic field h = 0.8 the expression 
in equation (17) has a zero at Ai « 0.788 + 0.337z, while the analytically continued function for the energy in the 
non-degenerate case has a branch cut which starts at a value of A equal to A2 ~ 0.898 + 0.217z and extends almost 
parallel to the x axis. Since y|Ai| — | A2 |, the branch cut of the square root in the degenerate case corresponds to 
the one of the non-degenerate case demonstrating the consistency of the two approaches starting from the classical 
ground state in the absence of a magnetic field and the magnetic field dependent classical ground state. 



IV. C20 



Next we consider the 20 site system on the vertices of a dodecahedron, shown in figure 19. In C20 there is only one 
kind of bond, which we call J, and each atom has three nearest neighbors as in the C12 case. In the classical solution 
spins belonging to the same pentagon are not coplanar jl3|], and the energy per bond is — = —0.7454. This is bigger 
than the energy per bond for the coplanar spins on an isolated pentagon, which is cos (^L) = -0.8090. Thus when 



the dodecahedron is assembled from the individual pentagons there is a cost in energy, in contrast to the twelve-site 
system case. This system has a discontinuity in < S z > as a function of the magnetic field at the classical level unlike 
the twelve-site system. The discontinuity has a magnitude of 0.64 and occurs at a magnetic field h = 1.432J. This 
is also the case for Cqq whereas in C12, C70 and Cs4 there is a discontinuity in the slope of the magnetization with 
applied field [ fl3"|| . We start again from the classical ground states and perturb them with the quantum fluctuations. 
The difference now is that the exact solution of the problem is not known and the Hilbert space is much larger with 
2 20 = 1,048,5 76 states. 

The question we wish to address is whether the differences between C12 and C20 seen at the classical level in 
magnetic properties survive for the S = \ case. As we will show below these differences do survive and are seen in 
the analytic structure of the perturbation theory for the two systems. Since spins are not coplanar at the classical 
level, the Hamiltonian, which is defined as described in section II, assumes the following complicated form: 

^ = H + \H 1 , (22) 



H Q — - — ^2 SfSj , (23) 

<ij> 

Hi= J2 faSfSt + P« S i S i + ^ S i~ S J + S v S ? S i) + h - c - > ( 24 ) 
<*,j> 

where a^, /3y, 7y and 6ij are complex coefficients defined analytically. For example, for the Si ■ S2 term a.\2 — 

j(l + -^), P12 = g(l — ir)(l ~ V^i), 712 = — g(l + V3i) and S12 = 7*3. Since here the starting point is the classical 
ground state, defining a local z axis at each site, the coefficients of Hi are in general different for different bonds. 
Furthermore they are complex because of the non-planar character of the spins in the classical ground state. We 
apply the same methods towards the solution as in the twelve-site system case. The elements of the 2x2 effective 
Hamiltonian matrix H e ^ in equation (10) are now real for all applied magnetic fields. 

A. Ground and Excited States 

The energies and wavefunctions of the two lowest lying states were calculated up to h — 0.72J. The ground state 
(o) and the excited state (o) are shown in figure 20. By drawing straight lines through the energies and extrapolating 
to zero field the lowest energy state in each spin sector is recovered, as well as the second lowest energy in the 5 = 
sector. We find these energies to be Ei{S = 0) = —9.722 J, E 2 {S = 0) = —9.345 J, E 3 (S = 1) = — 9.208J and 
E±{S = 2) = —8. 523 J. The magnetic field dependence is given by a Zeeman term —h^2i—i S'. Unlike C12 the first 
excited state at h < 0.137 is a singlet rather than a triplet |T3j. 

B. Correlation Functions 

In figure 21 we plot the lambda dependence of the expectation value < S z > for h = 0.4. We see a behavior similar 
to the one of C12, as the ground state assumes the value —2.7 x 10~ 8 and the excited state the value 0.9999997 at 
A = 1. Again this points to the accuracy of the calculated wavefunctions. The ferromagnetic or antiferromagnetic 
character of the correlation between any two sites remains the same in the ground state and the two lowest lying 
excited states at h — 0. In each of these states the nearest neighbor correlation ~ —0.3, while all other < Si ■ Sj > 
have magnitudes ~ 0.03 — > 0.08 at A — 1. Correlation functions are plotted for the lowest lying singlet and triplet 
states at h = 0.4, where they are the ground and first excited state in figures 22 and 23. We first consider the nearest 
neighbor correlations. For the singlet state, figure 22, the < Si ■ S2 > and < Si ■ S& > are equal only at A = and 
A = 1. The A dependence for < A < 1 depends on both the starting classical ground state and the value of h. For 
the lowest energy singlet state, < Si ■ S2 > and < Si • S§ > grow in magnitude from their classical value, —0.186 at 
A = 0, to —0.324 at A = 1 for all nearest neighbor correlations. Thus they are comparable to the value of the nearest 
neighbor correlation in the ground state of an isolated pentagon, —0.375. 

In the lowest lying S = 1 S z = 1 state, figure 23, two values of the nearest neighbor correlation functions are 
present, —0.297 and —0.327 at A = 1. This state can be described by considering the bonds on a pentagon, on its 
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diametrically opposite mirror and on a chain formed by the 10 remaining sites. On these bonds < Si ■ Sj >= —0.297 
while on the bonds connecting the chain to the two pentagons < Si ■ Sj >= —0.327. Since this configuration can be 
chosen in six different ways the triplet S — 1 S z — 1 state is six fold degenerate. There is little difference between 
the next nearest neighbors and more distant correlations in the singlet and triplet states. The quantum fluctuations 
significantly reduce the magnitude of these correlations compared with their classical values at A = 0. 

As a further check of the accuracy of the wavefunction the calculated values of these correlations give for < S 2 >= 
J2ij < &i ' Sj > th° value —0.00089 for the ground state and the value f. 99154 for the excited state at A = 1. These 
are consistent with the calculated < S z > values in figure 21. The change in the value of < S 2 > between the singlet 
and triplet states comes from changes in all the < Si ■ Sj >. The change due to the nearest neighbor < Si ■ Sj >'s 
is ~ 1 and the rest comes from small changes, ~ 0.005, from the other bonds. This suggests that there is no simple 
characterization of the excited states. 

In the classical ground state all the spins in the top and bottom pentagons have the same azimuthal angle a with 
the z axis, which is later taken to be the direction of the magnetic field. The rest of the spins form a different angle 
(3 with the z axis, so there are two different angles in the classical solution. There are only two distinct functions 
< Si ■ h > of the spins for the calculated states in the magnetic field, and they represent these two kinds of sites in 
the classical ground state. They were calculated for i = 1 and 6 and they are plotted in figure 24 for the ground and 
the excited state. They are consistent with < J2i Si ■ h >—< S z > h calculated in the ground and excited state in 
figure 21. 

For h < 0.137 the excited state is a singlet. This state has the same symmetry as in the triplet state with two 
types of nearest neighbor bonds equal to —0.304 and —0.324. As in the triplet case the longer range correlations, 
both ferromagnetic and antiferromagnetic are still comparable in magnitude and ~ 0.05. 

C. Analytic Structure 

Next we study the structure of the analytic continuation of the off-diagonal element B^ C (X). In figure 25 we plot 
the structure at a magnetic field equal to 0.3. Again we observe the asymmetry with respect to the y axis. As in 
the case of C12 the transition of the ground state between different spin sectors as a function of applied field can 
be tracked in the complex A plane. However the trajectory of the zeros of B% C (X), the off diagonal element in the 
effective Hamiltonian in equation (10), is different from C\i where the zeros move from the complex plane onto the 
real axis at A = 1 and h = h c and then move away for h > h c as seen in figure 17. These values of A are branch 
points of the energy function for the ground and first excited states. In C20 on the other hand close to the critical 
field, h c = 0.514, two zeros of B approach the real A axis from above and below at A ~ 0.85. The critical field is the 
value of h were the triplet state first becomes the ground state at A = 1. As h is increased further the two zeros move 
along the real axis in opposite directions, as shown in figure 26. At these values of A there are discontinuous jumps in 
the value of < S z > in the ground state as opposed to the rapid variation with A as seen in C\i . These jumps occur 
as the ground and excited states with different values of < S z > switch roles when B^ c - (A) changes sign. 

In figure 27 we plot < S z > for the ground state as a function of applied magnetic field for different values of A. 
The curves for A < 1 are strongly reminiscent of the results for the magnetization of the classical approximations 
for C20 and Ceo Ml - For ft, > 0.5 there are discontinuities in < S z > for different values of A. In these calculations 
starting from the field independent classical ground state the magnetic field is part of the perturbation and is scaled 
by A. So the effective field at which discontinuities in < S z > take place is Xh. In the classical approximation for 
C20 there is a jump in the magnetization as a function of h at 1.43 equal to 0.68S C , where S c is the magnitude of the 
classical spin. In the ground state calculated here the jump in < S z > occurs for A > 0.72 and grows in magnitude 
until, at A = 1, the jump discontinuity in < S z > is 1 for h c = 0.514. 

Different magnetic properties for the S — h solutions of C12 and C20 are thus seen to come from the different 
analytic behavior of the wavefunctions in the complex A plane. Similarity of the quantum solutions with the classical 
solutions suggests that this difference in analytic structure comes from the different connectivities of the two systems. 

V. CONCLUSION 

We have applied perturbation theory in the strength of quantum fluctuations around the classical ground states 
and found essentially exact results for the ground and first excited states for two frustrated spin systems, C12 and 
C2o- For C20 large orders in perturbation expansion and a high degree of numerical precision are required to get 



convergence of the analytic continuation of the matrix elements of the effective Hamiltonian and the wavefunction 
coefficients calculated perturbatively. 

We found that although sites in both systems are three-fold coordinated the spectra are qualitatively different with 
the lowest lying excited state being a triplet in C12 while it is a singlet in 620- The difference in the behavior of the 
two systems is attributed to their difference in the connectivity. The average energies per nearest neighbor bond for 
S = 5 spins, — 0.325J (in C20) and — 0.3105J (in C12), taking J\ = J2 = J, are comparable so that the "average" 
degree of frustration is the same for both systems. In C\i the first excited state arises by taking a linear combination 
of states each associated with a triplet excitation of a single non-triangle bond. By comparing < S\ ■ S4 > in figures 
14 and 15 the non-triangle bond is seen to take on a more strongly triplet character in the excited state. In C20 on 
the other hand there is no nearest neighbor bond which can be singled out to describe the character of the singlet and 
triplet excited states. The different values of S 2 for states in C20 arise from small changes in long range correlations. 
In order to distinguish between C12 and C20 we observe that in C20 all bonds are members of closed loops in which all 
bonds are equivalent while in Cyi the non-triangle bonds form alternating sides with triangle bonds to make hexagons. 
However characterization of connectivity by the nature of bonds forming closed loops is not sufficient in itself since 
the excited state of the Heisenberg model on a cube has a triplet excited state. 

In the classical approximation it was possible to characterize the ground state for different systems in terms of a 
topological number, the Skyrmion number, and to associate this with the presence or absence of discontinuities in the 
magnetization. We have shown that this difference, seen in the classical approximation for Cyi and C20, survives in S z 
for the exact states of the isotropic S = \ Heisenberg model and that it is associated with the analytic behavior in the 
complex coupling constant plane. Although it seems clear that it is the combination of frustration and connectivity 
which is responsible for the nature of the spectrum and the response of the system to a magnetic field, how exactly 
these are determined is not obvious. 

The connection between frustration and the character of the excitation spectrum is also of interest for infinite 



lattices. The apparent singlet nature of the excitation spectrum in the Kagome lattice arises from frustration [|18 26 
and leads to a large non magnetic contribution to the entropy of the system at low temperatures. This is an important 
feature in analysis of experimental data on any physical realization of a Kagome lattice. This result for the Kagome 
lattice may not generalize to other frustrated lattices if the sensitivity to connectivity seen in C12 and C20 is a guide 
to larger systems p7| . 

The accuracy of the wavefunctions found at A = 1 suggest that the Pade approximants derived from the perturbation 
series are very accurate representations of the wavefunctions and energies in the complex plane between the origin 
and A = 1 and that the analytic structure of these functions accurately represents the dependence on A and applied 
field. We have shown how the rapid change of < S z > for A ~ 1 and h ~ h c and the discontinuity in < S z > for C20 
can be traced back to the Pade approximants. The ability to track this order-disorder transition from coherent states 
to eigenstates of S z and S 2 as a function of A and applied field may be developed for richer many body systems with 
the calculational resources available today. 

These calculations are being extended to higher values of spin and to finite temperatures and will be used to 
investigate the temperature evolution of the properties of magnetic molecules, M1112, Fes etc., and, in particular, the 
quantum mechanical description of the large moments presently used in their analysis []28| . 
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FIG. 1. Space configuration of C12. Intra-triangle bonds are called Ji, while inter-triangle bonds are called J2. 




FIG. 2. Zeros and poles of the Pade approximant of the diagonal element An of the effective Hamiltonian for a = 1 in C12 
with 220 orders of perturbation theory used, o : zeros , x : poles. 
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FIG. 3. Zeros and poles of the Pade approximant of the diagonal element An of the effective Hamiltonian for a = 2 in C12 
with 220 orders of perturbation theory used, o : zeros , x : poles. 
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FIG. 4. Correlation functions for J2 = 1 in the ground state of C12 at h = : solid line : < Si ■ S2 > , dotted line : 
< Si • 54 >, dashed line : < Si ■ S5 > , long dashed line : < Si ■ Ss > 
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FIG. 5. Correlation functions for J2 = 1 in the first excited state of C12 at ft = 
< Si ■ 54 >, dashed line : < Si ■ §5 > , long dashed line : < Si ■ Ss > 
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FIG. 6. Nearest neighbor correlation functions for J2 = 2 in the ground state of C12 at h = : solid line : < Si ■ S2 > , 
dotted line : < Si ■ S4 >, dashed line : < Si • S5 > , long dashed line : < Si ■ Sn > 
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FIG. 7. Nearest neighbor correlation functions for J2 = 2 in the S — 1 , S z = first excited state of C12 at h = : solid line 
: < Si ■ S2 > , dotted line : < Si • S4 >, dashed line : < Si ■ S 5 > , long dashed line : < Si • S 8 > 
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FIG. 8. Ground state and first excited state energy as a function of applied magnetic field for C12, o : S z = 0, □ : S z = 1, 
o : S z = 2. 




FIG. 9. S* z vs. A for C12 for a magnetic field h ( ground state ) starting from the magnetic field independent classical ground 
state , J2 = 1.0 , solid-plus line : h=0.500 , dotted line : h=0.670 , dashed line : h=0.687 , long dashed line : h=0.688 , 
dot-dashed line : h=0.690 , solid line : h=0.700. 
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FIG. 10. S z vs. A for C12 for a magnetic field h ( excited state ) starting from the magnetic field independent classical 
ground state , J2 = 1.0 , solid-plus line : h=0.500 , dotted line : h=0.670 , dashed line : h=0.687 , long dashed line : h=0.688 
, dot-dashed line : h=0.690 , solid line : h=0.700. 




FIG. 11. S z vs. A for C12 for a magnetic field h starting from the magnetic field independent classical ground state , ,h = 1.0 
, solid line : h=1.010 , ground state , dotted line : h=1.010 , excited state , dashed line : h=1.050 , ground state , long dashed 
line : h=1.050 , excited state. 



1 n 




1.0 




0.0 0.5 1.0 

X 

FIG. 13. S z vs. A for C12 for a magnetic field h starting from the magnetic field independent classical ground state , J2 = 1.8 
, solid line : h=1.3423 , dotted line : h=1.400 , dashed line : h=1.500. 
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FIG. 14. Nearest neighbor correlation functions for C12 for J2 = 1 , h — 0.7 , ground state : solid line : < Si • S2 > , dotted 
line : < Si ■ §4 >, dashed line : < Si ■ S& > , long dashed line : < Si ■ S$ > 
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FIG. 15. Nearest neighbor correlation functions for C12 for J2 = 1 , h - 
line : < Si ■ S4 >, dashed line : < Si ■ §5 > , long dashed line : < Si ■ Ss 
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FIG. 16. Zeros and poles of the Pade approximant of the square root in C12 with 220 orders of perturbation theory used, 
magnetic field h=0.688, o : zeros , x : poles. 
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FIG. 17. Zeros of the Pade approximant for the square root using 270 orders of perturbation theory for various magnetic 
fields , o : h=0.686 , □ : h=0.687 , o : h=0.6875 , + : h=0.688 , x : h=0.689 , * : h=0.690. 
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FIG. 19. Space configuration of C20- All the bonds are equal to J. 
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FIG. 20. Ground state and first excited state energy for different values of applied magnetic field for C20 , : ground state 
, o : excited state. The dotted line is a fit for the ground state, and the long dashed line for the excited state. 
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FIG. 21. < S z > as a function of A for the ground state and first excited state of C20 for a magnetic field h = 0.4 starting 
from the magnetic field independent classical ground state , solid line : ground state , dotted line : excited state. 
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FIG. 22. Spin correlations in the ground state of C20 as a function of A for h=0.4: solid line : < Si ■ S2 >, dotted line : 
< Si- §6 >, dashed line : < Si ■ S3 >, long dashed line : < Si • Sg >, dot-dashed line: < Si • Sn >, thick solid line: < Si ■ Sis >• 
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FIG. 23. Spin correlations in the first excited state of C20 as a function of A for h=0.4: solid line : < Si ■ S2 >, dotted line : 
< Si • Se >, dashed line : < Si ■ S3 >, long dashed line : < Si • Sg >, dot-dashed line: < Si • Sn >, thick solid line: < Si • Sis >. 




FIG. 24. Ground and excited state spin correlations for h = 0.4 as a function of A for C20 at sites i=l and i=6: solid line : 
< Si ■ h > in the ground state, dotted line : < S 6 ■ h > in the ground state, dashed line : < Si ■ h > in the excited state, long 
dashed line : < S& ■ h > in the excited state. 




FIG. 25. Zeros (o) and poles (x) of the Pade approximant B^ C (X) in the complex A plane for C20 in an applied field h — 0.3, 
240 orders of perturbation theory were used. 
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FIG. 26. Zeros of the Pade approximant B^ c (\) along the A real axis for C20. 240 orders of perturbation theory were used 
and o and o are the values calculated. The dotted and long-dashed lines are fits. 




FIG. 27. < S z > as a function of magnetic field for the ground state of C20 for various A's : o : A = 0.70 , □ : A = 0.75 , + 
: A = 0.80 , x : A = 0.85. The lines are fits. 



